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Abstract 



A well-defined regularized path integral for Lorentzian quantum gravity in three 
and four dimensions is constructed, given in terms of a sum over dynamically tri- 
angulated causal space-times. Each Lorentzian geometry and its associated action 
have a unique Wick rotation to the Euclidean sector. All space-time histories pos- 
sess a distinguished notion of a discrete proper time. For finite lattice volume, the 
associated transfer matrix is self-adjoint and bounded. The reflection positivity of 
the model ensures the existence of a well-defined Hamiltonian. The degenerate geo- 
metric phases found previously in dynamically triangulated Euclidean gravity are 
not present. The phase structure of the new Lorentzian quantum gravity model can 
be readily investigated by both analytic and numerical methods. 
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In spite of numerous efforts, we still have not found a consistent theory of quan- 
tum gravity which would enable us to describe local quantum phenomena in the 
presence of strong gravitational fields. Moreover, the entanglement of technical 
problems with more fundamental issues concerning the structure of such a theory 
often makes it difficult to pinpoint why any particular quantization program has not 
been successful. There seems to be a broad consensus that a correct non-perturbative 
treatment should involve in an essential way the full "space of all metrics" (as op- 
posed to linearized perturbations around flat space) and the diffeomorphism group, 
i.e. the invariance group of the classical gravitational action, 



with k~ x = 87rGNewton and the cosmological constant A. 

One approach that does not rely on the existence of supersymmetry tries to 
define the theory by means of a non-perturbative path integral. The aim is not 
to evaluate this in a stationary-phase approximation, but as a genuine "sum over 
all geometries". Since even the pure gravity theory in spite of its large invariance 
group possesses local field degrees of freedom, such a sum must be regularized to 
make it meaningful. The two most popular approaches, quantum Regge calculus and 
dynamical triangulations 0, both employ simplicial discretizations of space-time, 
on which then the behaviour of metric and matter fields is studied. One drawback of 
these (mainly numerical) investigations is that so far they have been conducted only 
for Euclidean space-time metrics gf^f 1 instead of the physical, Lorentzian metrics. 
This is motivated by the analogy with non-perturbative Euclidean (lattice) field 
theories on a fixed, flat background, whose results can under suitable conditions 
be "Wick-rotated" to their Minkowskian counterparts. The amplitudes exp iS are 
substituted in the Euclidean path integral by the real weights exp(— S eud ) for each 
configuration. Real weight factors are mandatory for the applicability of standard 
Monte Carlo techniques and, more generally, for the convergence of the state sums. 

Unfortunately, it is not at all clear how to relate path integrals over Euclidean 
geometries to those over Lorentzian ones. This has to do with the fact that in 
a generally covariant field theory like gravity the metric is a dynamical quantity, 
and not part of the fixed background structure. For general metric configurations, 
there is no preferred notion of "time" , and hence no obvious prescription of how to 
Wick-rotate. 

One might worry that with a discretization of space-time the diffeomorphism 
invariance of the continuum theory is irretrievably lost. However, the example of 
two-dimensional Euclidean gravity theories provides evidence that this is not neces- 
sarily so. There, one can choose a conformal gauge in the continuum formulation, 
and obtain an effective action by a Faddeev- Popov construction. Physical quantities 
computed in this approach coincide with those computed using an intermediate dis- 
cretization in the form of so-called dynamical triangulations. In this latter approach 
one approximates the sum over all metrics modulo diffeomorphisms by a sum over all 




(1) 
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possible equilateral triangulations of a given topological manifold. Local geometric 
degrees of freedom are given by the geodesic lengths (all equal to some constant a) 
of the triangle edges and by deficit angles around vertices. Since different triangula- 
tions correspond to different geometries, the set-up has no residual gauge invariance, 
and it is straightforward to define (regularized versions of) diffeomorphism-invariant 
correlation functions. In the continuum limit, as the diffeomorphism-invariant cutoff 
a is taken to zero, the sum over triangulations gives rise to an effective measure on 
the space of geometries (whose explicit functional form is not known). 

This makes the regularization by dynamical triangulations, which is also appli- 
cable in higher dimensions, an appealing method for investigating quantum gravity 
theories. It is a further advantage that the formalism is amenable to numerical sim- 
ulations, which have been conducted extensively in dimensions 2, 3 and 4. Alas, in 
d = 3, 4, and for Euclidean signature, an interesting continuum limit has not been 
found. This seems to be related to the dominance of degenerate geometries. At 
small bare gravitational coupling GNewton, the dominant configurations are branched 
polymers (or "stacked spheres") with Hausdorff dimension <i# = 2, whereas at large 
^Newton the geometries "condense" around one or several singular links or vertices 
with a very high coordination number, resulting in a very large da- These extreme 
geometric phases and the first-order phase transition separating them are qualita- 
tively well described by mean field calculations |2|, [J. Another unsatisfactory aspect 
of the Euclidean model (as well as other discretized models of gravity) is our inability 
to rotate back to Lorentzian space-time. 

In order to tackle these problems, two of us have recently constructed a Lorentz- 
ian version of the dynamically triangulated gravitational path integral in two dimen- 
sions [|J. The basic building blocks are triangles with one space-like and two time-like 
edges. The individual Lorentzian geometries are glued together from such triangles 
in a way that satisfies certain causality requirements. The model turns out to be 
exactly soluble and its associated continuum theory lies in a new universality class 
of 2d gravity models distinct from the usual Euclidean Liouville gravity. One central 
lesson from the two-dimensional example is that the causality conditions imposed 
on the Lorentzian model act as a "regulator" for the geometry. Most importantly, 
they suppress changes in the spatial topology, that is, branching of baby universes 
is not allowed. As a result, the effective quantum geometry in the Lorentzian case 
is smoother and in some senses better behaved: a) in spite of large fluctuations 
of the geometry, its Hausdorff dimension has the canonical value g?# = 2, unlike in 
Euclidean gravity, which has a fractal dimension du = 4; b) in spite of a strong 
interaction between matter and gravity when the system is coupled to Ising spins, 
the combined system remains consistent even beyond the c—1 barrier, unlike what 
happens in Euclidean gravity [|], |]. 

Motivated by these results, we have constructed a discretized Lorentzian path 
integral for gravity in 3 and 4 space-time dimensions. Unlike in two dimensions, the 
action is no longer trivial, and the Wick rotation problem must be solved. Also the 
geometries themselves are more involved, and the geometry of the (d— l)-dimensional 
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spatial slices is no longer described by just a single variable. We have succeeded in 
constructing a model with the following properties: 

(a) Lorentzian space-time geometries (histories) are obtained by causally gluing 
sets of Lorentzian building blocks, i.e. <i-dimensional simplices with simple 
length assignments; 

(b) all histories have a preferred discrete notion of proper time t; t counts the 
number of evolution steps of a transfer matrix between adjacent spatial slices, 
the latter given by (d— l)-dimensional triangulations of equilateral Euclidean 
simplices; 

(c) for a fixed space-time volume Na (= number of simplices of top-dimension), 
both the Euclidean and the Lorentzian discretized gravity actions are bounded 
from above and below; 

(d) the number of possible triangulations is exponentially bounded as a function 
of the space-time lattice volume; 

(e) each Lorentzian discrete geometry can be "Wick-rotated" to a Euclidean one, 
defined on the same (topological) triangulation; 

(f) at the level of the discretized action, the 'Wick rotation" is achieved by an an- 
alytical continuation in the dimensionless ratio cn = —l 2 ime /l 2 pace of the squared 
time- and space-like link length; for a — — 1 one obtains the usual Euclidean 
action of dynamically triangulated gravity; 

(g) for finite lattice volume, the model is (site) reflection-positive, and the transfer 
matrix is symmetric and bounded, ensuring the existence of a well-defined 
Hamiltonian operator; 

(h) the extreme phases of degenerate geometries found in the Euclidean models 
cannot be realized in the Lorentzian case. 

For the sake of definiteness and simplicity, we will concentrate mostly on the 
three-dimensional case. The discussion carries over virtually unchanged to d = 4, the 
details of which will be given elsewhere JTJ. (Obviously, the corresponding continuum 
theories will be very different, one describing a topological quantum field theory, and 
the other a field theory of interacting gravitons.) The classical continuum action 
is simply eq. (|l|), with d = 3. Each discrete Lorentzian space-time will be given 
by a sequence of two-dimensional compact spatial slices of fixed topology, which 
for simplicity we take to be that of a two-sphere. Each slice carries an integer 
"time" label t, so that the space-time topology is I x S 2 , with I some real interval. 
The metric data will be encoded by triangulating this underlying space by three- 
dimensional simplices with definite edge length assignments. There are two types 
of edges: "space-like" ones (of length squared I 2 = a 2 > 0, with the lattice spacing 
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a > 0) which are entirely contained in a slice t — const., and "time-like" ones (of 
length squared I 2 = —aa 2 < 0) which start at some slice t and end at the next slice 
t + 1. This implies that all lattice vertices are located at integer t. 

A metric space-time is built up by "filling in" for all times the three-dimensional 
sandwiches between t and t + 1. We only consider regular gluings which lead to 
simplicial manifolds. Our basic building blocks are given by three types of Lorentzian 
tetrahedra, 

(1) 3-to-l tetrahedra (three vertices contained in slice t and one vertex in slice 
t + they have three space- and three time- like edges; their number in the 
sandwich [t, t+1) will be denoted by N 3 i(t); 

(2) l-to-3 tetrahedra: the same as above, but upside-down; the tip of the tetra- 
hedron is at t and its base lies in the slice t+1; notation iy 13 (t); 

(3) 2-to-2 tetrahedra: one edge (and therefore two vertices) at each t and t + 1; 
they have two space- and four time-like edges; notation N 22 {t). 

Each of these triangulated space-times carries a discrete causal structure ob- 
tained by giving each time-like link an orientation in the positive t-direction. Then, 
two lattice vertices are causally related if there is a sequence of positively oriented 
links connecting the two. 

The discretized form of the Lorentzian action ([I]) is obtained from Regge's pre- 
scription for simplicial manifolds f|, see [|7| for details. The action is written as a 
function of the deficit angles around the one- dimensional edges and of the three- 
dimensional volumes, which in turn can be expressed as functions of the squared 
edge lengths of the fundamental building blocks. The contribution to the action 
from a single sandwich [t, t+1] is 

AS a (t) = 4nak y ^+(N 31 (t)+N 13 (t))(akK 1 ~a 3 \L 1 ) + N 22 (t)(akK 2 -a 3 \L 2 ), (2) 

where 

K x {a) 
K 2 (a) 
H<*) 

Note that the sandwich action (|j) already contains appropriate boundary contribu- 
tions, such that S is additive under the gluing of contiguous slices. In relation (0), 
we have used the rescaled cosmo logical constant, A = kA. 

At each time t the physical states \g) are parametrized by piece- wise linear ge- 
ometries, given by unlabelled triangulations g of S 2 in terms of equilateral Euclidean 
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V3a+T V2a + 1 



triangles. For a finite spatial volume N (counting the triangles in a spatial slice), the 
number of states is exponentially bounded as a function of N and the vectors \g), 
defined to be orthogonal, span a finite-dimensional Hilbert space TCn- The transfer 
matrix T/v will act on the Hilbert space 

N 

H (N) := Hi, (4) 

?'= TV ■ 

1 * mm 

and the states \g) will be normalized according to 

(gM = -^S gi>92 , ^2C g \g){g\ = l. (5) 

U 9i g 

The symmetry factor C g is the order of the automorphism group of the two-dimen- 
sional triangulation g, which for large triangulations is almost always equal to 1. 
With each step At = 1 we can now associate a transfer matrix T/v describing the 
evolution of the system from t to t+1, with matrix elements 

(g2\f N (a)\ gi ) = G a (g u g 2 ; 1) = £ e iAS «. (6) 

sandwich (gi —>g2 ) 

The sum is taken over all distinct interpolating three-dimensional triangulations of 
the "sandwich" with boundary geometries g\ and g2, each with a spatial volume 
< N. The propagator Gj^(gx, g 2 ;t) for arbitrary time intervals t is obtained by 
iterating the transfer matrix t times, 

G N (g 1 ,g 2 ;t) = (g 2 \f t N \g 1 ), (7) 

and satisfies the semigroup property 

G N (g 1 ,g 2 ; h + t 2 ) = J2 C 9 G N (g2, g\ h) G N (g, gr, h). (8) 

9 

The infinite- volume limit is obtained by letting iV — > oo in eqs. 

A brief remark is in order on our notion of "time" : the label t is to be thought of 
as the discretized analogue of proper time, as experienced by an idealized collection 
of freely falling observers. We do not claim that this is a distinguished notion of time 
for pure quantum gravity, but it is a possible choice, in the present case suggested 
by our regularization. Note that in continuum formulations the proper time gauge 
is not usually considered, because it is a gauge choice that - considered for arbitrary 
geometries - goes bad in an arbitrarily short time. This problem does not occur in 
the discrete case: by construction we only sum over space-time geometries for which 
there is a globally well-defined (discrete) "proper time" . 

The action S associated with an entire space-time S 1 x S 2 of length t in time- 
direction is obtained by summing expression (fj) over all t' = 1, 2, . . . t and identifying 
the two boundaries. The result can be expressed as a function of three geometric 



6 



"bulk" variables, for example, the total number of vertices Nq, the total number of 
tetrahedra A3 and t, 



S a (N , A 3 , t) = A (AakiK! - K 2 ) - 4a 3 A(Li - L 2 )) + A 3 (akK 2 - a 3 AL 2 ) 

+t Uak(-K^fa - 2(Kx - K 2 )) + 8a 3 A(Li - L 2 )) . (9) 



Because of the well-known inequality iVo < (A 3 + 10)/3, valid for all closed three- 
dimensional simplicial manifolds, this implies the boundedness of the discretized 
Lorentzian action at fixed three-volume. This result is analogous to what happens 
in Euclidean dynamical triangulations. We write the partition function as 

Z a (k,\,t)= e lS ^ T ^ T ^\ (10) 

TeT t (S 1 xS 2 ) 

with %{S l x S 2 ) denoting the set of all Lorentzian triangulations of S 1 x S 2 of length 
t. It will turn out that a necessary condition for the existence of a meaningful contin- 
uum limit is the exponential boundedness of the number of possible triangulations 
as a function of the space-time volume A 3 : only if the growth is at most expo- 
nential in A 3 , can this divergence potentially be counterbalanced by a cosmological 
constant term exponentially damped in A 3 . In our case, exponential boundedness 
follows trivially from the same property for Euclidean triangulations (where it has 
been proven rigorously for d = 3,A || since the Lorentzian space-times form a 



subset of the former. Note that the convergence of the partition function implies 
the absence of divergent "conformal modes". 

As it stands, the sum (|10|) over complex amplitudes has little chance of converg- 
ing, due to the contributions of an infinite number of triangulations with arbitrarily 
large volume A3. In order to make it well-defined, one must perform a Wick rotation, 
just as in ordinary quantum field theory. Thanks to the presence of a distinguished 
global time variable in our model, we can associate a unique Euclidean triangu- 
lated space-time with every Lorentzian history contributing in (|S|). It is obtained by 
taking the same topological triangulation and changing the squared lengths of all 
time-like edges from — aa 2 (Lorentzian) to aa 2 (Euclidean), leaving the space-like 
edges unchanged. We can then use Regge's prescription for calculating the (real) 
Euclidean action S®(N , N 3 , t) associated with the resulting Euclidean metric space- 
time (where a is always taken to be positive). After some algebra one verifies that 
by a suitable analytic continuation in the complex a-plane from positive to negative 
real a, the Euclidean and Lorentzian actions are related by 

S_ a (N ,N 3 ,t)=iS*(N ,N 3 ,t), (11) 

for a > \ . This latter inequality has its origin in the triangle inequality for Euclidean 
triangles: S E is real only for a > |, and the limit ct = \ corresponds to the degenerate 
case of totally collapsed triangles. Moreover, a = — 1 is the only point on the real axis 
in which the coefficient of t in the Lorentzian action (Op vanishes (this corresponds 
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to a — 1 in eq. (JTTJ) ) . In this case one rederives the familiar expression employed in 
equilateral Euclidean dynamical triangulations, namely, 



= Sf = -ak(2irNi - 6N 3 arccos -) + a 3 AiV 3 ^. (12) 

Our strategy for evaluating the partition function is now clear: for any choice of 
a > |, continue (|]) to —a, so that 

e iS a (N ,N 3 ,t) <x-*-a e ~ S% (N ,N 3 ,t) _ ^g-j 

Because of the exponential boundedness, the Wick-rotated Euclidean state sum in 
([Hp is now convergent for suitable choices of the bare couplings k and A. We can 
therefore proceed in two ways: either attempt to perform the sum analytically, by 
solving the combinatorics of possible causal gluings of the tetrahedral building blocks 
(as has been done in d = 2 Q), or use Monte-Carlo methods to simulate the system 
at finite volume. Once the continuum limit has been performed, we can rotate back 
to Lorentzian signature by an analytic continuation of the continuum proper time T 
(which in the case of canonical scaling is of the form T—at; not to be confused with 
the triangulation T) to iT . If we are only interested in vacuum expectation values 
of time-independent observables and the properties of the Hamiltonian, we do not 
need to perform the Wick rotation explicitly, just as in usual Euclidean quantum 
field theory. 

Let us now establish some properties of the discrete real transfer matrix T = 
T(a=—1) of our model that are necessary for the existence of a well-defined Hamil- 
tonian, defined as h := — ^T 2 . These will be useful in any proof of the existence 
of a self-adjoint continuum Hamiltonian H. It is difficult to imagine boundedness 
and positivity arising in the limit from regularized models that do not have these 
properties. 

We will show that Tjy is symmetric, bounded for finite spatial volume N, and 
that the two-step transfer matrix T 2 is positive. Symmetry is obvious by inspection. 
The sandwich action (§) is symmetric under the exchange of in- and out-states 
(corresponding to AT 31 <-> AT 13 ). So is (j^), because the counting of interpolating 
states does not depend on which of the geometries g±, g% defines the in-state, say. The 
boundedness of T/v for finite spatial volume N follows from the finite dimensionality 
of the Hilbert space Ti.^ it acts on and the fact that there is only a finite number of 
possibilities to interpolate between two given spatial triangulations g\ and g<i in one 
step. Positivity of the two-step transfer matrix, > 0, follows from the reflection 
positivity of our model under reflection with respect to planes of constant integer-t 
(for regular lattices, this property is also referred to as site-reflection positivity, 
c.f. |P-1]|). In order to be able to define a Hamiltonian as 

h N ■= -j^f 2 , (14) 
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we must make sure that the square of the transfer matrix is strictly positive, > 0. 
We do not expect that T N has any zero-eigenvectors, because this would imply the 
existence of a "hidden" symmetry of the discretized theory. It may of course happen 
that there are "accidental" zero-eigenvectors for certain values of N. In this case, 
we will adopt the standard procedure of removing the subspace J\f( N > spanned by 
these vectors from the Hilbert space, resulting in a physical Hilbert space given by 
the quotient H$P = H^/Af^K 

It should be emphasized that although the summation in the path integral is 
performed in the "Euclidean sector" of the theory, our construction is not a priori 
related to any path integral for Euclidean gravity proper. The point, already made in 
the two-dimensional case ||, is that we sum only over a selected class of geometries, 
which are equipped with a causal structure. Such a restriction incorporates the 
Lorentzian nature of gravity and has no analogue in Euclidean gravity. We therefore 
expect our Lorentzian statistical mechanics model to have a totally different phase 
structure from that of Euclidean dynamical triangulations. This expectation is 
corroborated by an analysis of the "extreme phases" of Lorentzian quantum gravity, 
to determine which configurations dominate the path integral 

Z*(k,\,t)= £ e"^, (15) 

TeT t (S 1 xS 2 ) 

for either very small or very large inverse Newton's constant k > 0. In order to 
make a direct comparison with the Euclidean analysis [|T2| , [|, we set without loss of 
generality a — 1 in eq. (15) and rewrite the Euclidean action (12) as 



Sf = k 3 N 3 -k 1 N 1 , (16) 



with the couplings 



ki=2iiak, k 3 = 6 ak arccos — \-a 3 X — =. (17) 

In the thermodynamic limit N 3 — > oo, and assuming a scaling behaviour such that 
t/N 3 — > 0, one derives kinematical bounds on the ratio of links and tetrahedra, 
f := Ni/N 3 , namely, 

1<£<^- (18) 

This is to be contrasted with the analogous result in the Euclidean case, where 1 < 
£ < |. It implies that the branched-polymer (or "stacked-sphere" ) configurations, 
which are precisely characterized by £ = |, and which dominate the Euclidean state 
sum at large k±, cannot be realized in the Lorentzian setting. The opposite extreme, 
at small k±, is associated with the saturation of the inequality 

iVi < In (N - 1), (19) 
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and in the Euclidean theory goes by the name of "crumpled phase". At equality, 
every vertex is connected to every other vertex, corresponding to a manifold with a 
very large Hausdorff dimension. Again, it is impossible to come anywhere near this 
phase in the continuum limit of the Lorentzian model. Instead of fll9D , we have now 
separate relations for the numbers and iVj of space- and time-like edges, 

N[ sl) = £(3iVo(i) - 6) = 3N - 6t, Jvf } < £ N (t)N (t + 1). (20) 
t t 

Assuming canonical scaling, the right-hand side of ineq. (|19"D behaves like (length) 6 , 
whereas the second relation in ( ETj ) scales only like (length) 5 . 

We conclude that the phase structure of Lorentzian gravity must be very differ- 
ent from that of the Euclidean theory. In particular, the extreme branched-polymer 
and crumpled configurations are not allowed in the Lorentzian theory. This is an- 
other example of causal structure acting as a "regulator" of geometry. It also raises 
the hope that the mechanism governing the phase transition will be different, and 
potentially lead to a non-trivial continuum theory, in three as well as in four dimen- 
sions. 
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